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ABSTRACT 



One of the important light elements created during the big bang nucleosyn- 
thesis is 7 Be which then decays to 7 Li by electron capture when recombination 
becomes effective but well before the Saha equilibrium recombination is reached. 
This means that 7 Be should wait until its recombination epoch even though 
the half-life of the hydrogenic beryllium atom is only 106.4 days. We calculate 
when the conversion from primordial 7 Be to 7 Li occurs taking into account the 
population of the hyperfine structure sublevels and solving the kinetic equations 
for recombination, photoionization and conversion rate. We also calculate the 
energies and the spectrum of narrow neutrino doublet lines resulting from 7 Be 
decay. 
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Introduction 



The first major event i n the Universe after the big bang and the inflation stage is 



the creation of light elements (I Wagoner. Fowler, and Hoyldll967t ) including D, He-3 and Li 



which are of importance for our everyday life or considered as energy source on the Earth in 
the future. The combination of neutron and proton to form deuterium is the first step in the 
nuclear reaction chain. The big bang nucleosynthesis (BBN) starts when the temperature of 
the Universe is about T ~ O.lMeV (redshift z ~ 4 x 10 8 ) and the deuterium destroying high 
energy photons have become rare. BBN proceeds swiftly after this deuterium bottleneck is 
overcome and is over by T ~ 0.01 Me V. The next major event in the Unive r se is tha t of th e 
recombination of primordial plasma (IZeldovich. Kurt, and Syunyaevi (119681 ) jPeebled (119681 )) 



to form neutral atoms starting with helium III to helium II recombination at 4500 < z < 7000 
and ending with hydrogen recombination at 800 < z < 1800. In this paper we will explore 
a unique connection between these two major events, the recombination of beryllium 7 Be 
and its decay to lithium 7 Li and in the process also mention a doublet of narrow lines in 
the cosmic neutrino background. This unique event brings together the MeV scale nuclear 
physics of nucleosynthesis and eV scale atomic physics of recombination. We will also find 
that the very low energy atomic physics of hyperfine transitions in 7 Be has a surprisingly 
strong influence on the decay properties of 7 Be usually associated with MeV scale physics. 
This problem of beryllium to lithium conversion is one of the last problems not investigated in 
detail so far connected with the primordial nucleosynthesis and the process of recombination 
of ions and electrons in the radiation field of the cosmic microwave background (CMB). This 
problem is very attractive because it demonstrates a very rare case in astrophysics when 
the atomic processes and even the populations of the hyperfine structure influences nuclear 
processes. 

The cosmological parameters used for numerical calculations are Hubble constant 
Hq = 73km/s/Mpc, matter fraction Q m = 0.24, CMB temperature Tcmb = 2.72QK, baryon 
fraction = 0.043, cosmological constant Q\ = 0.76. 
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2. Decay of 7 Be 



At the end o f BBN the abundances of light elements relative to hydrogen are 
f lSerpico et al\\20<m : X 4He = 6 x 10" 2 , X 2H = 2.5 x 10~ 5 , X 3He = 10" 5 , Xa H ~ 1(T 7 , 
Xr Bc ~ 10~ 10 , Xr Li ~ 10 -11 , Xe u ~ 10~ 14 and X M < 10~ 15 , where X M is the abundance 
of all heavier elements, Xi = rij/nn, is the abundance of species i and tir- is the number 
density of hydrogen nuclei. Note that about ~ 90% of the primordial 7 Li comes from 7 Be 
and only ~ 10% of 7 Li is p roduced directly in BBN, this division is however sensitive to the 
baryon to photon ratio rj (llocco et a/.ll2009l ). 



7 Be in its fully ionized state in the early Universe is stable. Fully ionized b eryllium 



can c apture electrons from the plasma at the high densities in the stellar interiors (jBahcall 



19891 ). However the electron capture rate for the fully ionized beryllium from the plasma is 
completely negligible in the early Universe due to the low electron density. Once recombined 
with an electron beryllium decays by electron capture through the following reactions 



7 Be + e~ -»■ 7 Li + u e 
7 Be + e - -> 7 Li* + z/ c , 7 Li* ->■ 7 Li + 7, 



(1) 
(2) 



where 7 Li* is the excited state of lithium nucleus with energy 477. 6KeV above the ground 
state. The Q- value of the reaction is 861.8KeV (ITilley et a/.ll2002l ). The laboratory value 
for the half life of 7 Be is 53.2 days. The branching ratio is 10.44% for reaction [2J These 
laboratory values have been calculated and measured for neutral 7 Be with four electrons. 

The 7 Be in the primordial plasma will however decay as soon as it acquires a single 
electron. Therefore we are interested in hydrogen-like 7 Be 3+ . The electrons in the 2s shell 
of neutral 7 Be have a ~ 1% effect on the decay ra te as show n by experiments involving 
Be 2+ (OH 2 )4 as well as other chemical compositions (jHuhlll999l ) and we will ignore the in- 
fluence of 2s shell electrons in the discussions below. Helium-like 7 Be 2+ has practically the 
same lifetime and branching ratio as the neutral 7 Be. The change in the decay rate and the 
branching ratio from helium-like 7 Be 2+ to hydrogen-like 7 Be 3+ depends strongly on the spin 
temperature of the hyperfine structure of hydrogen-like 7 Be 3+ ions. 
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;o understand. We will make simple estimates of decay rates 



This statement is easy 

based on the calculations by IPatyk et al\ (]2008[ ). 7 Li* has a nuclear spin of I = 1/2 which 
is different from the nuclear spin of / = 3/2 of 7 Li and 7 Be. The two initial states available 
for 7 Be 3+ are with total angular momentum F + = 2 and F~ = 1. The two final states 
available to 7 Li* + z/ e are F = 0,1. Thus the angular momentum can only be conserved 
for the reaction [2] from the initial state F~ = 1 of 7 Be 3+ to the final state F — 1 of 
7 Li*. The reaction is suppressed for the other initial hyperfine state with total angular 
momentum F + = 2. Correspondingly for reaction (TJ where the nuclear spin does not change 
in the reaction, the two available final states have the same total angular momentum as 
the initial state, F = 1,2, and the reactions from both the initial hyperfine states are 
possible. Millielectronvolt hyperfine splitting dictates the branching ratio for these two 
nuclear reactions. In hydrogen-like 7 Be 3+ ion the population levels of the two hyperfine 
states (or equivalently the spin temperature) are determined by the rates of radiative decay 
and pumping due to well known physical processes. 

Our calculation is simplified because pumping by the CMB radiation field and due to 
electron collisions is much faster than the decay rate and cosmological time, i.e. the spin 
temperature should be close to the CMB temperature T equal at that time to the electron 
temperature. The energy of the hyperfine splitting is much smaller than the temperature 
at that time. This means that the spin temperature is very high compared to the hyperfine 
splitting and the hyperfine structure sublevels are populated according to their statistical 
weights. 



The Lyman-a resonant scattering (IField 



1958 



Varshalovich 



19671 ) with rate P a 



47r J dvo v l v j (his) IPOs -1 at z ~ 30000, spin chan ging collisions with electrons with rate 
~ n e /«io ~ 10~ 4 s _1 ( Furlanetto and Furlanettoll2007 ) and stimulated emission and absorption 
of CMB at hyperfine transition with rate ~ Bi\lv ~ 4 x 10 _6 s _1 are all faster than the decay 
rate A ~ 7 x 10 _8 s~ 1 and would make the spin temperature equal to the matter/radiation 
temperature T. Above a v is the absorption cross section for Lyman— a photons including the 
line profile from Doppler broadening. I v is the background CMB intensity. It corresponds to 
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Wien region of the black body spectrum for the Lyman-a photons and Rayleigh- Jeans region 
for the hyperfine transition radiation. B21 is the Einstein B coefficient for stimulated emis- 
sion. n e is the el ectron number density and Km is the s pin change cross section. n w oc Z~ 2 for 
hydrogenic ions (jAugustin. Miiller. and Greinerl ll989). where Z is the nuclear charge. There 
is also a mild energy dependence of the cross section at temperatures of interest of ~ 1/T 1 / 2 . 
The spin cha nge cross sect i on for hydrogen-electron co l lision s is of the order ~ 10~ 9 cm 3 s _1 
at T ~ 10 5 K (ISmithl Jl966f ) ]Furlanetto and Furlanettol (120071 )) and thus for 7 Be 3+ it will be 
~ 10 _10 cm 3 s _1 . Thus Lyman-a scattering, Ionization and recombination, Collisions with 
electrons and stimulated emission and radiation of CMB and maintains T spin = T. These 
rates are plotted in Figure [1] along with the expansion rate defined by the Hubble parameter 
H(z). It is interesting to note that a typical hydrogenic beryllium atom would have flipped 
its electron spin a billion times due to Lyman-a scattering and recombined and ionized 100 
times before finally being able to decay. The energy difference between the two hyperfine 
states is 10 _4 eV << T. Thus the hyperfine states would always be distributed according to 
their statistical weights. 

In helium-like 7 Be 2+ (or Li-like and neutral 7 Be) spin directions of two Is shell electrons 
is always opposite and the problem radically differs from the problem in the hydrogen-like 
ion. Thus there is no hyperfine splitting in the case of helium-like ions, but the relative 
direction of the spins of the captured electron and nucleus is of great importance. We can 
write the probabilities of the two reactions in helium-like ions in terms of the corresponding 



probabilities o: 
calculations of 



the hydrogen-like ions in the two hyperfine states using the recent theoretical 



Patyk et al. 



( 12008) which have the supp ort of experiments involving hydrogen- 
like and helium-like 142 Pm ions (IWinckler et a/.ll2009l ). Then we can use these relations and 
the relative populations of the hyperfine structure sub-levels in the hydrogen-like ions to find 
the formula for the probabilities of both channels of interest. The reactions for hydrogen-like 
7 Be 3+ and helium-like 7 Be 2+ are summarized in Figures [2] and |3j 



For reaction [T] the initial and final spin of the nuclei are same, I = 3/2 and the electron 
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Fig. 1. — Comparison of different processes influencing beryllium to lithium decay. It is 
interesting to note that a typical hydrogenic beryllium atom would have flipped its electron 
spin a billion times due to Lyman-a scattering and recombined and ionized 100 times before 
finally being able to decay. 



capture rates of 7 Be 2+ and 7 Be 3+ are related by (IPatyk et a/.ll2008l ) 

'2F+ + 1\ ^ (IF- + 1 



A 



7Be 2 + 



A + + 



(3) 



21 + 1 J V 2J + 1 
where X* 1 are the decay rates from two hyperfine states of 7 Be 3+ . If f ± are the fraction of 
atoms in the two hyperfine states, the net decay rate is given by 



Ar B e3 + = / + A + + f-X 

( 2F+ + 1 

= v 2(27+1; 

= ^7 Bc 2 + /2, 



A+ + 



2F~ 



2(2/ + 1) 



(4) 



where we used the fact that (2F + + 1) + (2F~ + 1) = 2(2/ + 1). The factor of two can be 
understood as follows. In the helium-like atom both channels corresponding to A + and A - 
are available simultaneously which are added according to the statistical weights. For the 
hydrogen-like atom also for high spin temperature we add the probabilities of two channels 
according to the statistical weights. But for the helium-like atom the capture rate gets 



7 R 3 + F=1 
Be 

F=2 



Q=861.8KeV, T 1/2 = 106.4 days 

Fig. 2. — Effect of hyperfine splitting of 7 Be 3+ on the nuclear reactions. The half-life given 
above corresponds to the high redshift universe when the spin temperature is maintained at 
a high value by collisions with electrons. The half-life would be different if the populations of 
hyperfine levels of hydrogenic 7 Be 3+ are not distributed according to their statistical weights. 

multiplied by a factor of two corresponding to the fact that any of the available electrons 
can be captured through any of the available channels. Thus the factor of two is the result 
of the fact that there are two electrons available for helium-like atom and that for high 
spin temperature for the hydrogen-like atom we add the capture rates according to their 
statistical weights. This factor would, in general, be different from two if the atoms in the 
hyperfine state were not distributed according to their statistical weights, for example at 
low temperatures comparable to the energy difference between the two states. Since the 
statistically averaged rate coefficients relevant for our calculation for the helium-like and 
hydrogen-like atoms differ by a just factor of two, we do not need the individual rates for 
A + and A~ and can use the experimental result for the helium-like/neutral 7 Be. 

For reaction [2] however the 7 Li* has a spin of 1/2 and only capture from F~ state is 
allowed. The above formula is valid with A + = and for the net decay rate we have the 
same result, A* Be3+ = A* Be2+ /2. We are using * for the reaction rates of reaction [2] while no * 
indicates the reaction rates for reaction[TJ Using the laboratory value of Ai a b = ln2/53.2days 
for the total decay rate of 7 Be 2+ , we get net decay rate Abbn = 0.5Ai a b with the same branching 
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Li + v e 

Q=861.8KeV, T 1/2 =53.2 days 

Fig. 3. — Decay of helium-like 7 Be 2+ . 



ratio, 10.4% of the decays following reaction [2J The half life of hydrogen-like 7 Be 3+ with the 
hyperfine states distributed according to their statistical weights is thus equal to 106.4 days. 
Coincidently this is the same result we would have got from the simple considerations of the 
electron density at the nucleus. We must emphasize here that the exact life time of 7 Be 3+ 
is not very important for us as long as it is much shorter than the cosmological time at the 
epoch of beryllium to lithium conversion and the duration of 7 Be 3+ recombination. All the 
width of the two neutrino lines comes from the kinetics of recombination. In particular we 
have ignored the electron-electron interaction in the helium-like 7 Be 2+ , the effect of which 
on our calculation of time of 7 Be decay and the neutrino spectrum would be ~ few%. 



3. Recombination of 7 Be 

Now we can write down the kinetic equations for recombination of 7 Be 4+ . 

[n e (z)X7 Be 4+a7 Be 4+ - /37 Be 3+X7 Be 3+] (5) 



dXTQ e i+ 1 



dz H(z)(l + z) 

■ [— n e (z)X 7Be i+a7 Be i+ + /37 Be 3+X7 Be 3+ + A b bn^7 Be 3+] (6) 



dX7 Be 3+ 1 



dz H(z)(l + z) 
X 7Be = Xrl^ + X 7Be 4+ + X 7Be 3+ , (7) 
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where n e (z) is the number density of electrons at redshift z, H(z) is the Hubble parameter, 
«7 Be 4+ is the total recombination coefficient including recombination to the ground state. 
j^decay j g ^ e lithium fraction coming from decay of beryllium and excludes the lithium pro- 
duced during BBN. Because of the extremely low number density of beryllium, the number 
of ionizing photons released during direct recombination to the ground state is negligible 
(about 10 orders of magnitude less) compared to the ionizing photons already present in the 
background radiation and can be neglected. /37B e 3 + is the total ionization coefficient which is 
related to «7 Be 4+ by the condition that Saha equation must be satisfied in equilibrium. Thus 

{27im e k B Tfl 2 

/3 7Be 3+ = a7Be 4+ e l bt . (8) 



Note that the 2-photon decay rate from 2s state of 7 Be 3+ is 3.4 x lOV 1 flGoldmaJl9891 ). This 
is much faster than the same rate for hydrogen ~ 8s _1 due to the higher charge of the nucleus. 
More importantly it is much faster than the ionization rate ~ 6 x 10 _3 s _1 at z ~ 30000 
and the electron capture rate X^ n = 7.5 x 10~ 8 s -1 , £i/ 2 = 106.4days. Thus 2s level must 
be included in the total recombination coefficient. The number of Lya photons produced 
would be of order X7 Be and can be neglected. Thus equations for beryllium recombination 
are much simpler than that of hydrogen and helium recombination. «7 Be 4+ is the total 
rec ombination coefficient including direct re combinations to the ground state and is given 
by (IPequignot. Petitjean. and Boissonlll99ll ) 

a 7Be , + = cmV\ (9) 

where t = (T/10*K)/Z 2 , a = 5.596, b = -0.6038, c = 0.3436 and d = 0.4479 and Z is 
the nuclear charge.. Figure |4] shows the result of integrating the recombination equations 
for beryllium. For comparison the Saha equilibrium solution for beryllium recombination is 
also shown. The beryllium to lithium conversion occurs significantly earlier at z — 30000 
than the z = 25000 value predicted by the Saha solution. The reason of the difference is 
connected with the short decay time of recombined 7 Be 3+ . In the Saha equation we follow the 
balance between the recombination and photoionization but a typical atom has recombined 
and ionized many times even though the net recombination in equilibrium may be small. In 
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reality due to the decay of beryllium on a time scale much shorter than the cosmological 
time the equilibrium is never established. 
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Fig. 4. — Lithium number density as a fraction of total beryllium and lithium number density 
from decay of beryllium. For comparison thermal equilibrium results using Saha equation 
are also plotted. 



4. Energy release and narrow doublet of neutrino lines 



Every decay of beryllium to lithium is accompanied by emission of an electron neutrino 
with 89.6% of them having energy of 861.8KeV and nuclear recoil energy of 56.5eV. 10.4% of 
the decays go to excited state of lithium, electron neutrino gets 384.2KeV and nuclear recoil 
has energy of only 11.2eV which is slightly larger than the thermal energy T = 7eV of the 
plasma at redshift z = 30000. The recoil velocity of v/c = 5.85 x 10 -5 will result in a natural 
line width of 56eV due to Doppler shift for the 477.6KeV photon emitted when lithium 
relaxes to ground state almost instantaneously (7i/ 2 = 73fs). These photons will down 
scatter by C ompton scattering with electro ns giving most of their energy to the plasma due to 
recoil effect ( Zeldovich and Sunyaev fl969 ). The energy transfer by Compton scattering will 
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become inefficient when the photons reach the critical energy of m P H (z)/re P crT L=annnn ~ 80eV 



when the energy transfer rate becomes less than the expansion rate ( iBernstein and Dodelson 



( 119901 ) ) . This will leave a small distortion in the high energy part of CMB. The energy 
transfered to plasma will additionally cause a y— type distortion of ~ (0.104i? 7 /T)X7Be^ ~ 
10~ 16 . The neutrinos however will free stream to us and we can calculate the neutrino 
spectrum today. This is plotted in Figure [5] including both 861.8KeV neutrinos and 384.2KeV 
neutrinos, 



0.896(1 + zf dX^ ay 0.104(1 + zf dX%* y 

1 i „ 861.8KeV H 

861.8KeV dz ll+2= b 384.2KeV dz 



(10) 

where njjo is the density of hydrogen nuclei at z = 0. The full width at half maximum 
(FWHM) for the first line is 2.3eV and the central energy is 29.5eV. For the second line 
the FWHM is leV with a central energy of 13.1eV. The line width at half maximum is 
SE/E = 7.8%. The width and asymmetric line profile of these lines is defined by the 
kinetics of recombination. 



5. Comparison with other sources of neutrinos and energy release before 

recombination 

At T ~ O.lMeV annihilation of electrons and positrons makes the biggest contribution 
to cosmic neutrinos, however these neutrinos have a much broader thermal spectrum and 
form part of the low energy cosmic neutrino background (CNB). Decay of neutrons and other 
nuclear reactions during BBN also contribute neutrinos, but due to high redshift z ~ 10 8 
these neutrinos are also today redshifted to the sub-eV energies, ~ lMeV/10 8 ~ O.OleV. The 
spectrum of these neutrinos would be broad since they are emitted over a period when redshift 
changes by a factor of ~ 2. Energy released during BBN is also quickly thermalized and does 
not lead to spectral distortions. Decay of tritium created during BBN occurs at z ~ 2.5 x 10 5 
corresponding to half life of 12.32 years. Decay of tritium to helium-3 results in an electron 
with average energy of Q e = 5.7KeV and an antineutrino with average energy ~ 12.9KeV. 
The antineutrino spectrum would be broad as they decay over a period corresponding to 
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Fig. 5. — Neutrino spectrum from beryllium decay in the early Universe. There are two lines 
corresponding to two electron capture decay branches. For comparison equilibrium results 
(marked Saha) are also plotted. The FWHM for high energy line is 2.3eV with a central 
frequency of 29.5eV. For second line FWHM is leV central frequency is 13.1eV. The line 
width at half maximum is 5E/E = 7.8%. The width and asymmetric line profile of these 
lines is defined by the kinetics of recombination. 

the age of the Universe at that time and these antineutrinos would have an average energy 
of ~ 0.05eV today assuming they have zero mass. About 10% of tritium has decayed by 
z = 6.35 x 10 5 and 90% by z=1.35 x 10 5 leading to the width of the neutrino spectrum of 
SE/E ~ 1.4. This broadening is comparable to the intrinsic width of the neutrino spectrum 
SE/E ~ E max / E avg ~ 18.6/12.9 = 1.4. Thus non-thermal neutrinos from all sources before 
recombination other than 7 Be decay would have a spectrum much broader than those from 



7 



Be decay and would have much lower energy. These 



relati vistic today in more than one mass eigenstates (INakamura and Particle Data Group 



ow energy neutrinos would be non- 



20101 ) compared to relativistic neutrinos from 7 Be decay. The energy released into the plasma 
from tritium decay in the form of energetic electrons would result in a chemical potential 
of \i ~ (Q e /T)(n3^/n 7 ) ~ 10~ 15 in the CMB. This is about the same amount of entropy 
generated in the beryllium decay much later. 
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Conclusions 



Lithium- 7 observed today was original 



nucleosynthesis (I Wagoner. Fowler, and Hoyle 



y pro duced as beryllium-7 during primordial 



19671 ). Although half life of beryllium atoms 



is very short, it has to wait until the beginning of recombination epoch to decay. We have 
calculated the exact redshift when this happens. We have also estimated the effect of energy 
release during the decay on the cosmic microwave background. In addition the neutrinos 
produced during the decay give rise to unique narrow lines in the cosmic neutrino background. 
These lines are too weak to be observable today. We should mention that the detec tion of the 



much more numerous but lower energy CNB neutrinos is currently being discussed ( jRingwald 



20091 ). The recombination and decay of beryllium has nevertheless theoretical significance. 



It marks the end of primordial nucleosynthesis and the beginning of recombination. 
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